Spatially resolved spectroscopy on a superconducting proximity nanostructure 
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We investigated the local density of states (LDOS) of a normal metal (N) in good electrical contact 
with a superconductor (S) as a function of the distance x to the NS interface. The sample consists 
of a pattern of alternate stripes of Au and Nb made by UV-lithography. We used a low temperature 
scanning tunneling microscope to record simultaneously dl /dV{V, x) curves and the topographic profile 
z{x). Nearby the NS interface, all the spectra show a dip near the Fermi energy but depending on the 
geometry, different behaviours can be distinguished. First, when th e characteristic size of the normal 
metal L is much larger than the coherence length ^jv — y/hDN the spectral extension of the dip 
decreases from A at the NS interface to zero at distances a; ^ ^jv. Second, when L is comparable to 
^jv the apparent gap in the LDOS is space-independent and related to the Thouless energy. 



PACS numbers: 74.50.+r 73.40. Gk 73.50.Bk 
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A normal metal in good metallic contact with a super- 
conductor can acquire some superconducting properties 
and reciprocally the superconductivity can be affected by 
the normal metal vicinity. This phenomenon is known as 
the proximity effect. It has been first extensively stud- 
ied in the late 1960's using the Ginzburg-Landau the- 
ory which describes macroscopic superconducting sys- 
tems near their transition temperature Tc Recent 
technical and material developments renewed the interest 
in the physics of this effect, especially at the mesoscopic 
scale. Simultaneously a more comprehensive understand- 
ing of the proximity effect in the diffusive regime based 
on the theory of non-equilibrium superconductivity has 
emerged For instance, predictions on the spatial de- 
pendence of the local density of states (LDOS) of a prox- 
imity structure are made. At very low temperature, when 
the characteristic size L of the normal metal becomes 
smaller than the thermal length Lt = ^TiDn /2'n:kBT 
(_Djv is the diffusion constant in the normal metal, ks 
the Boltzman constant, A the superconducting gap), the 
variation of the LDOS of a NS structure is predicted to 
depend on the ratio L /^jv with ^at = \/'hDj^ the co- 
herence length When L ^ £^n, the superconducting 
correlations extend to distances larger than leading 
to a depression of the electronic density of states around 
Ep-, the Fermi energy. Experimental tests of this theory 
were obtained by Gueron et al. 0. They used nanofab- 
ricated tunnel junctions to probe the electronic density 
of states n{E, x) at three distinct positions along a nor- 
mal wire in contact with a superconductor. Scanning 
tunneling microscopy (STM) was also used to study this 
proximity effects in ballistic N metals {i.e. L <^ S^j^ 
where Z^r is the elastic mean free path in N) §,|. In 



these experiments the induced gap crucially depends on 
the thickness of the normal system. The shape of the 
spectra is explained within the de Gennes-Saint James 
bound states model |jl^. In a diffusive sample where 
Zat <C i < ^AT, a gap whose value depends on the Thou- 
less energy, Eth = hD/L^, is predicted STM 
experiments on small N wires embedded in a S matrix 
showed a healing length of superconductivity on the S 
side much larger than the coherence length ^g. But ac- 
cording to the authors, they could not easily interpret the 
spatial dependence of the LDOS in the N metal because 
of a complicated geometry. 

In this Letter we report measurements of the local elec- 
tronic density of states (LDOS) by scanning tunneling 
spectroscopy on Nb/Au proximity junctions. We spa- 
tially resolve the LDOS in the normal metal as a function 
of the distance x to the NS interface. We can discrimi- 
nate between the two situations: L < and L » 
within the same sample. This rules out any spurious 
effects due to modifications of the film properties and 
interface qualities. 

Spectroscopy by STM allows a high energetic and spa- 
tial resolution in conjonction with sample morphology. 
The STM hangs inside a sealed tube by a one meter 
long spring to decouple it from external vibrations. This 
tube is immersed in a He* cryostat. Cooling is achieved 
by introducing high purity helium exchange gas into the 
tube. The temperature is reduced to 1.5 K by pumping 
on the helium bath. Tunneling spectra are obtained with 
a lock-in detection technique with a 77 /iV peak-to-peak 
modulation voltage of the bias at IkHz. We measure 
the differential conductance dl/dV versus V while hold- 
ing the STM tip at a fixed height above each position 
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FIG.y (a) Optical microscope picture of the sample. It con- 
sists of a 50 nm thick pattern of self aligned stripes 1 nm wide 
and 250 fim long of Au and Nb. The line represents section of 
the film. Inset (b) is a scanning electron microscope photog- 
raphy of the sample displaying the small ridges on each side 
of the Au stripe. Part (c) is a set of two spectra taken above 
either Au or Nb. The experimental curves have been shifted 
for clarity. Above the Nb film (lower curve) the spectrum is 
fitted within a BCS model with A = 1.15 meV and T = 0.07A 
(solid line). In contrast spectra taken above Au (upper curve) 
can not be fitted within this model. The fit cannot describe 
simultaneously the wings and the zero-bias conductance. The 
dotted line is a fit with A — 1.1 meV, F = 0.12 A whereas the 
solid one is for A = 1.1 meV, F = 0.40A 

r. This provides a local probe of n{eV,r). The energetic 
resolution is determined by the thermal broadening at 
1.5 K. Each curve is a single acquisition process with no 
additional averaging. The dl/dV curves and the topog- 
raphy are measured during the same line scan allowing 
us to correlate them precisely with the distance x to the 
NS interface. Data are normalized to the conductance 
at high voltage {V > A). The same normalization co- 
efficient has been used for all the curves. The typical 
tunneling resistance is 10^ fl. 

We have chosen Au as a normal metal because it is 
chemically inert. The Au film has been previously char- 
acterized by measuring the temperature dependence of 
the resistivity, R{T), which gives a mean free path In ~ 
22 nm and a diffusion coefficient Dn ~ 1.0 x 10~^ m^s~^ 
leading to ^jv — 53 nm. Nb is used as a superconductor 
because of its high critical temperature. A 50 nm thick 
Nb ffim undergoes a BCS transition at 8.1 K and shows an 
energy gap A ~ 1.15 me V at 1.5 K. Its coherence length 
is £,s = \/^Ds/2A ~ 27 nm and the mean free path is 
Is — 6nm given by previous R(T) measurements. The 
ratio Zs/fs < I is characteristic of a dirty superconduc- 
tor. Our sample consists of a 50 nm thick pattern of 
juxtaposed stripes, 1 /xm wide and 250 /im long, of Au 
and Nb in good electrical contact, see fig. 0(a) and|l|(b). 

The Nb was first uniformly deposited by DC sputter- 
ing. We used UV-lithography and subsequent reactive 
ion etching (RIE) to produce the Nb stripes. We then 
sputtered 2nm Ti and 50 nm Au and lifted-off the re- 
maining resist. As seen in figure0(c), Au forms ridges 



on each side of the Nb stripes. The size of these ridges 
has been measured both by STM and scanning electron 
microscopy. Their lateral dimensions vary from 50 nm 
to 200 nm and their height from 10 nm to 50 nm. They 
are in electrical contact with Nb but poorly connected 
to the rest of Au. Indeed we observe BCS-like dI/dV{V) 
curves in between the Au ridges and the Au stripes which 
indicates the presence of bare Nb (see lower curve of 
fig|l|(c)). The solid line represents the convolution of the 
density of states with a Fermi-Dirac distribution at the 
measured temperature T = 1.5 K. We introduce a phe- 
nomenological parameter F to define a complex energy 
E* — E — iV [|l^. This parameter reflects the broad- 
ening of the energy levels due to the finite lifetime of 
the quasi-particles. The calculation provides a super- 
conducting gap A = 1.15 meV and a finite lifetime of 
the quasi-particles F = 0.07A. In contrast spectra taken 
above Au either over the ridges or the stripes can not 
be fitted using this model. As an illustration, the upper 
part of fig. 0(c) shows several curves with different sets 
of parameters. The transition from these two types of 
spectra allows us to infer the NS interface. 

In order to describe our results we need to take into 
account coherence effects of the quasi-particles inside the 
normal metal. In the framework of the quasi-classical 
Green's function formalism, induced correlations between 
electrons of opposite spin in the normal metal are de- 
scribed by a complex function e{E, r) The LDOS 
of the quasi-particles n is related to 6 by n{E, r) — 
nf)Re[cos9{E, r)]. In the dirty limit and in zero magnetic 
field, 9{E, r) obeys the ID-diffusion Usadel equation [ p5[ : 

^ 1^ + [{~iE + F,„) - 2F,/ cos 9] sin 9 

-hA(a;)cos6' = 0, (1) 

where T sf and Tin are the spin-flip and the inelastic 
scattering rates respectively, A(a;) the pair potential. In 
the normal metal [x > 0) we assume no pair potential 
(Atv = 0) whereas in the superconductor {x < 0), A{x) 
obeys a self-consistent equation. For simplicity, we take 
A equals to its bulk value everywhere in the supercon- 
ductor. Conditions of continuity, 9^=0- = 9^=0+ and 

spectral current conservation (Ts(ff)£c=o- = <yN{^)x=o+ 
are imposed at the NS interface in the case of a perfect 
electrical contact. 

When L ^ ^7v, the normal metal can be considered as 
semi-infinite and one recovers the N density of states far 
away from the superconductor. The boundary condition 
of eq.(|l|) on the N side is then 9{E, oo) ^ 9n ^ 0. The 
resulting LDOS n{E,x) is predicted to display a space- 
dependent pseudo-gap structure Q. This situation cor- 
responds to the 1 fim wide stripe of Au. In the first set 
of experiments we scanned above its surface and mea- 
sured the LDOS as a function of the distance x to the 
NS interface. The spectra were taken along the same 
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FIG.|: (a) 3D plot of dl/dV versus V curves for different po- 
sitions of tlie tip above tlie wide stripe of Au. Tlie curves are 
evenly shifted along the x axis (not to scale) . The pseudo-gap 
decreases from A at the NS interface to at distances larger 
than 200 nm. (b) represents the normalized zero-bias conduc- 
tance (open circles) and the energy of the peak (full squares), 
Smax, as a function of the distance to the NS interface. The 
solid lines correspond to the fit. On the right hand side, we 
compare measured spectra (c) with those calculated (d) for 
three distances: 11 nm, 68 nm and 185 nm. 
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FIG.^: (a) shows a 400nmx400nm STM topographic image 
of a ridge 20 nm high. (b)is a measured line profile, (c) the 
3D plots represents spectra for different positions above the 
Au ridge, (d) the zero bias conductance, Gv=o/Gv=4,m.v , and 
the energy of the peaks, Smax are plotted as a function of the 
distance x to the NS interface. Smax ~ 1-5 meV does not vary 
while Gv^o/Gv=4,mV continuously increases with x. 



line perpendicular to the interface. The space dependent 
LDOS is plotted in fig.H(a) for different x. 

We observe a peak in Gv{x) — dI/dV{V, x) at an en- 
ergy Smax which decreases as we move away from the 
superconductor. Figure^(b) displays Smax and the nor- 
malized zero bias conductance Gv=o/Gv=4mV as a func- 
tion of X. A zero density of states at the Fermi level is 
never observed due to pair breaking mechanisms such 
as spin-flip and inelastic scattering which cannot be ne- 
glected at 1.5 K. We compare our measurements to the 
convolution of the LDOS calculated from eq. (|l|) with a 
Fermi-Dirac distribution. For the superconducting part, 
we used the value of A and Tin given by the BCS fit 
above Nb. For the normal part, the rates Tin, ^sf and 
the ratios $,s/£,n, <^s/'^n are taken as adjustable param- 
eters. We find Tin = 0.01 A, Tgf = 0.15A and the ratios 
S,s/S.N — 0.2, as/o'N = 4. The deduced mismatch pa- 
rameter 7 = (T'N^s/o'sS.N equals 0.05 instead of 2.5 eval- 
uated from the conductance measurements obtained on 
films on the same thicknesses. This discrepancy reflects a 
much stronger induced superconductivity in the normal 
metal than expected ||l^. This effect may be due to in- 
termixing between Au and Nb at the interface or to the 
presence of Ti in between these two materials. Indeed 
the transport measurements do not provide any informa- 
tion on the influence of Ti on the conductivity of Au. 
Moreover from Tin and Tsf we can deduce a decoherence 
time T — h/max{Tin, Tgf) ~ 20 ps. This value is in good 



agreement with recent experiments measuring the phase 
relaxation time on Au films deposited above either Ti 
0orAlll|].^ 

The fluctuations of Go and Smax are not due to experi- 
mental noise but are correlated to changes in the Au film 
morphology related to different grains. A slight disper- 
sion of the properties of these grains may explain such 
non-monotic evolutions 

In the second step of experiments, we have investi- 
gated the LDOS over two small ridges of Au. They are 
lying over Nb and the NS interface is now horizontal. 
Due to the initial slope of their edges, when moving the 
STM tip horizontally we can change the distance to the 
NS interface (see figures § and ^. Unlike previous ex- 
periments where the proximity effect was probed on NS 
bilayers with different thicknesses we do not change 
the overall size of the normal system between different 
spectra but only the distance to the NS interface. There- 
fore we obtain the actual spatial dependence of the LDOS 
n{E, x) of a normal system with a fixed dimension L. 

In this configuration, we find that the position of the 
peaks in dI/dV{V,x) does not depend on the distance 
to the NS interface in striking contrast to the preced- 
ing case. For a ridge of 20 nm height, we observe that 
the maximum is located at 1.5 meV (see fig.||) whereas 
in a ridge of height 50 nm, the position of the peak is 
0.95 meV (see fig.§). Theoretically, in such a finite dif- 
fusive geometry i.e. when L is of the order of ^at, the N 
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FIG.^: Picture (a) represents a STM topographic image of 
550 nmx550 nm, we have plotted a profile of the measured 
ridge. Spectra have been taken above this ridge along the 
white line in (a). The arrows represents places where spectra 
have been taken. For clarity dl /dV{v) curves taken from left 
to right have been shifted in (c). The ones correspponding 
to the black arrows are all similar except for the one corre- 
sponding to the white arrow taken in the bottom of the ridge. 
We have plotted Smax and Gv^o/Gv=4,mV for all the similar 
spectra as a function of x in (d) to prove that they are space 
independent. 

side boundary condition of eq. is ^{E,x = L) = 
and a minigap Eg is predicted to appear in the LDOS. 
Eg is related to the Thouless energy Eth = hDi^/L^ and 
Eg = mm(A, A(l-HL/2.1,^Ar)-2) ||^|. in the first ridge, 
Erh = 18meV^ A, the minigap is limited by the su- 
perconducting one. On the contrary, in the second ridge, 
Exh = 2.8 meV that yields to a minigap i5g=0.53meV. 
Eventhough, we measure the position of the peaks and 
not the value of the minigap, our results are consistent 
with the expected behaviour of the gap as a function of 
the size L. Further numerical simulations are still needed 
to get a more quantitative comparison. 

The evolution of the zero bias conductance with x is 
different for the two ridges. In the bigger one, the LDOS 
suddenly changes as soon as the tip probes the Au surface 
and remains identical for any location above the ridge 
(fig.^ whereas in the smaller one, Gv=QlGv=imV con- 
tinuously increases with x from 0.1 to 0.9 (fig. This 
can be due to 3D geometrical effects since in this lat- 
ter case the normal part is in a quasi-ballistic regime 
{L~In). 

In conclusion, we have probed the proximity effect in 
two different systems. First, we have investigated the 
LDOS of a semi-infinite normal metal and found space 
dependent energy spectra as a function of the distance 
to the NS interface. This behaviour is in good agreement 
with the pseudo-gap model predicted by the theory of 
non-equilibrium superconductivity. We have also inves- 
tigated the local density of states in a confined geometry. 



In this case we have found a spectral structure which 
does not vary in space and which can be related to the 
Thouless energy. However the evolution of the zero bias 
conductance depends on the shape of the normal metal. 

We wish to acknowledge beneficial discussions with M. 
Sanquer and the valuable help of J.C. Toussaint for nu- 
merical calculations. 
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